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SUMMARY 



The laminar -boundary- layer thickness, the boundary-layer Reynolds 
number and minimum critical Reynolds number, and the roughness Reynolds 
number have been calculated by an approximate method for a sphere and 
disk in the supercritical Reynolds n umb er region. The calculations for 
the sphere show that the boundary layer at the stagnation point of a 
sphere Is much thicker than that on an airfoil, that the boundary-layer 
thickness increases very slowly with an increase in distance from the 
stagnation point, that the boundary layer over the forward portion of a 
sphere Is highly stable at large Reynolds numbers with respect to the 
To l 1ml en-Schlichting type of waves, and that roughness of a given height 
produces the largest roughness Reynolds numbers at about 57° from the 
stagnation point. The calculations for the disk show the unusual result 
that the boundary-layer thickness Is greatest at the stagnation point 
and decreases with an increase in distance from this point, that the 
boundary layer is extremely stable with respect to the Tollmien- 
Schlichting type of waves, and that roughness of a given height produces 
a given roughness Reynolds number over a smaller portion of the disk 
surface than over the sphere surface. 


INTRODUCTION 


In connection with an experimental subsonic-speed investigation 
undertaken recently at the Langley Aeronautical Laboratory of the transi- 
tion from laminar to turbulent flow on a sphere at a Reynolds number of 

4.7 X 10 6 , based on a radius of 60 inches, the thickness of the laminar 
boundary layer, its stability, and the roughness Reynolds numbers for 
various-sized roughness particles were calculated by approximate methods 
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for Incompressible flow. In addition to the calculations • for a sphere, 
calculations were also made for a much blunter body, namely, for a disk 
with a radius of 60 inches perpendicular to the flow at a Reynolds number 

of 4.7 x 10 6 . 

Some time after the calculations were completed it was discovered 
that measured transition positions on a 9-inch-radius sphere had been 

reported in reference 1 up to a Reynolds number of 1*75 X 10^, based on 
the radius. Calculated distributions of the boundary- layer momentum 
thickness over the forward portion of the sphere for six Reynolds numbers 

f. 

between 0.5 X 10° and 3 X 10 as well as the position of the neutral- 

stability point up to a Reynolds number of 3 X 10^ were also presented 
In reference 1. 

In the present paper the calculated distribution of the boundary- 
layer momentum thickness over the forward portion of a sphere Is pre- 
sented in the form of a single nondimens ional curve that a^Lows the 
thickness at any Reynolds number to be obtained with almost no calcula- 
tion. The calculated boundary-layer Reynolds numbers at any value of 
the sphere Reynolds number can also be readily obtained from a single 
curve. The distribution of the minimum critical Reynolds number on the 
sphere, which is independent of the sphere Reynolds number. Is presented 
as well as the location of the neutral- stability point for Reynolds num- 
bers up to 1,000 x 10^. Finally, the distribution over the sphere of a 
nondimensional roughness parameter that allows the roughness Reynolds 
number to be quickly calculated for any small roughness height at any 
value of the sphere Reynolds number is shown In the form of a single 
curve. 

The calculated results are for the supercritical Reynolds number 
range, that is, for the range in which the final separation of the bound- 
ary layer on the sphere is that of the turbulent boundary layer. In this 
range the pressure distribution over the forward portion of the sphere is 
almost Independent of the Reynolds number, and, thus, the present calcula- 
tions are good approximations for any sphere at any supercritical Reynolds 
number. 

Most of the calculated quantities presented for the sphere are also 
presented in the same form for the disk. 


SYMBOLS 


a,A,b,C,d,e constants 


c 


reference length, radius of sphere or disk 
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F' 


velocity ratio , 


u/U (eq.. ( 15 ) of ref. 2) 


t 8 

f = 

kLU 



c 


h height of roughness particle 

K Schlichting velocity-profile shape parameter 


v dx 


Z 

dx 


V 


-Lt R 

u e 





reference Reynolds number, 
roughness Reynolds number. 



e,c 


mln-tnmTn critical Reynolds number, value of Rg at which a 
small disturbance is neither damped nor amplified 



c 


r 


perpendicular distance from point on surface of body of 
revolution to axis of symmetry 



m 

U 


velocity at outer edge of boundary layer and in direction of 
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U free-stream velocity 

00 


u = 

u 


u 




velocity inside boundary layer and in direction of x 
velocity at y = £ with roughness particle absent 


w 


2V 

Ucc 


W 


velocity through surface, positive outward 


x = * (For a sphere, x/c is equal to the angle in radians from the 
c stagnation point) 



increment in x 
initial station of x 

distance along surface, measured from stagnation point 
nondimensional distance from surface (ref. 2) 

distance normal to surface, positive outward 


3 wedge-angle parameter 

?! = Zr 2 


8 displacement thickness, 8/c 

5 * “ I 0 ( x - IK 


measure of boundary-layer thickness (ref . 3) 
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0 momentum thickness, 5/c 

* = So i(" ' §) d? 

[I viscosity of fluid 

v kinematic viscosity, p/p 

p density of fluid 

t surface shearing stress 

Subscripts: 
a actual 

d disk 

initial point 
plate 

stagnation point 
theoretical 
v surface 

Barred quantities are dimensional. 

ANALYSIS AND DISCUSSION 

Sphere 

The boundary-layer characteristics for a sphere were calculated by 
use of the modified Schlichting method (ref. 4) which was altered to 
apply to a body of revolution. (See ref. 3 for the original method.) 
The change from the method for two-dimensional flow consists in using 
the momentum equation for a body of revolution instead of the momentum 
equation for two-dimensional flow. The momentum equation for the flow 
over a body of revolution (see ref. 3) is 


o 

P 

t 

th 
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d| + 5 dU L + O ! 8 dr ^w _ ^ 

dx U dx\ 0/ r dx U pU c 


When the following three quantities 

d 2 - - 

k = £ au = (if -i = e 2 ££ R 

v dx \cj ^ x v dx c 


z sa 


(z = e\) 


L = -V„8R . tJE 

1 V U„ C V . 


T 0 T 

■P - w _ W R 

ST.’^ 9 

from reference 4 are introduced into equation (l) , the result is 




Introducing the quantity = Zr 2 into equation ( 2 ) gives 


In order to find {^(x) , equation ( 3 ) can he written as 

h “ 2 X* T f - *1 - 4 - v) + p i,o W 

where S, =0 because B-. — Zr 2 t and at x = 0 the quantity Z is 

l,o -l- 


finite and r ■ 0. 
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In order to integrate equation (4), Simpson's rule (see ref. 6) was 
used to obtain 



Equation ( 5 ) can be solved by iteration j the first approximation for 

+ ^-) l is obtained by extrapolation of 

9/ 4lx 

p-^ from the values at smaller values of x. Near the stagnation point, 

stiff icient values for applying this method of integration are obtained by 
the process described in the immediately following paragraph. When p^ 

and, thus, Z are known at any value of x, the quantities f, kq, k, 

and 5*/ 0 at that value of x can be calculated or read from charts . 
(See refs. 3 and 4.) 

In order to begin the integration at the stagnation point, equa- 
tion (2) is written as 



and use is made of the fact that, near the stagnation point, r = x and 
U = 0^0 x. The term ^ in equation (6) then becomes 

or 


ZU dr = 
r dx 


zi 


- z s 


the expression 


-Jd[f - - k(: 
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Equation (6) can then be written as 


i - §[ f - k i - 4 + £) - *t] 



u 


which, upon noting that U = ( — ) x near x = 0, gives 

\dx/ t 


,/az\ 

14 = 


lim 

x->0 


4 - ^ - 4 + zl ~ 


u 


Because U = 0 at x = 0, it is necessary that 


[ f - ■'i - 4 + t )\ ■ 0 


in order that 


be finite. 


Equation (8) together with the relation. 


g 2 (K +1) - f - k = 0 


(7) 


(8) 


(9) 


from reference 4 can be used to determine the initial values of K, k, 

and k n when v . , ( 4^ ] , and R„ are known by the same method used 
1 w,t 7 \<3x/t c 

for two-dimensional flow in reference 4. The quantity K is the param- 
eter that determines the shape of the velocity profile) g is a ratio 
of two boundary-layer thicknesses and depends on . K only. In the present 
calculations k]_ = 0, and, thus, equations (8) and ( 9 ) reduce to 


£ 


S (K + 1) 



( 10 ) 


Because f , g, and 5/9 depend only on the velocity-profile shape 
parameter K, equation (10) can be solved numerically for Kj. and the 

value of k-j- can then be calculated from equation ( 9 ) with k^_ = 0. The 
results are 


Kj. = -0.72542 

kt = 0.05299 
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Hie initial values of K^. and kj. are, thus, known; for the present 
case = 0 and, therefore, either K^. or ^ is sufficient to deter- 
mine K, f, S*/0, and g. 


In order to begin the integration of equation ( 4 ), it is useful to 
know how j3, varies near x = 0 . This variation is obtained from equa- 

tion ( 3 ) together with r 
*( 0 are finite. 

at x = 0. The second derivative 


k, and 5 
dx 


= x, U = (^j^J fact that f. 

Then, it follows from equation ( 3 ) that 
d 2 £l 


* 1 ' 


at x = 0 can also 


dx £ 


be easily calculated; thus, near x = 0 equation (3) becomes 



Then, 
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P x = -2^. — [f - - k^2 + ^)J + (Error) o[x^] (ll) 

(dx)t 


•Hie calculation of Pp is begun by malting use of equation (ll) and is 

■z 

continued by use of equation (5) together with the relations U = — sin x 


and r = sin x. (See p. 417 of ref. 7 .) 


Pi 


The quantity Z is obtained from P-. by the relation Z = — 
which follows from the definition of P-^. By using the definition 
k = Z equation (ll), the relation Z = where r = x, and equa- 

CLX 

tion (8), the subsequent relation results at the stagnation point: 

*t = 7 (12) 



All the calculations were ma de for a velocity distribution outside the 

3 

boundary layer given by U = ^ sin x,. the potential-flow velocity dis- 
tribution for a sphere. A comparison of the velocity distribution given 
by the expression U = ^ sin x with the measured velocity distributions 
(ref. 8) for the so-called supercritical range of Reynolds numbers shows 
that the expression U = ^ sin x is a good approximation for the experi- 
mental velocity distribution for most of the forward portion of the sphere 
and can be considered as the limiting distribution for the experimental 
velocity distributions. 
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The calculated variation of (or /z) with x/c is shown 

c ' c 

in figure 1. The presentation of the results in the form of & plotted 


against x/c enables the boundary-layer momentum thickness to be easily 
calculated for a sphere of any size and at any Reynolds number in the 
supercritical Reynolds number region. It is noted that the boundary- 
layer thickness increases very slowly with increase in distance from the 
stagnation point. This behavior differs from that on an airfoil (see 
p. 211 of ref. 9) on which the boundary- layer thickness increases rapidly 
from its value at the stagnation point. It is also noted that the value 

of at the stagnation point of the sphere is about 10 times as 

large as that at the stagnation point of the airfoil having an NACA 
64-A010 section at zero angle of attack. 


In figure 2 is shown the variation with x/c at E c = 4.7 x 10^ 
of the boundary-layer Reynolds number Rg and of the boundary-layer 
minimum critical Reynolds number Rq c , the Reynolds number below which 

very small wave-like disturbances ( Tollmlen-Schlichting type of waves as 
discussed, in ref. 9) are damped out. The Reynolds number Rg Is obtained 

from the relation 


U 9 PmS 
Uoo c v 

= U9R C 

= (us)( e #r) to) 

Because U and 9 ^R c are independent of R c , it is apparent that 
Rg is directly proportional to \^R^. Therefore, the value of Rg at 
any value of R c Ib obtained by using 


Rg = / ^ 

^ Rq )r c =4.7xio 6 1(4-. 7 x 10^ 


to) 


The minimum critical Reynolds number R fi depends on the shape param- 

0 > c 


eter K alone of the velocity profile. The values 
fore, obtained from figure 1 of reference 4 


of Rc 


are , the re - 


8,c 

The parameter K is depend- 
ent only on the velocity distribution outside the boundary layer and is 
independent of R c ; therefore, Rg c is independent of R c . Here, again, 

is different from that on an airfoil. 


the behavior of 


l c » 


and R 


8,c 
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(See ref. 4.) For the sphere the value of Rg Increases linearly with 
x/c for the smaller values of x/c. The values of R Q remain almost 

' o,C 

constant at first and then decrease at an increasing rate. For the air- 
foil with an NACA 64A010 section and at zero angle of attack, R 0 
increases in a parabolic manner with xj and Rg c , after a sudden large 
drop from its value at the stagnation point, remains almost constant to 
the value of x where the pressure is at a minimum. 

From equation (14) and figure 2 the value of R c necessary for 
neutral stability ^Rg = Re,c) an y value of x/c can be calculated. 

The results are shown in figure 3 • It is noted that at a Reynolds num- 
ber of h .7 X 10^ the boundary layer is stable to about x/c * 1.25, or 72° 
from the stagnation point. In order to make the boundary layer unstable 
as far forward as 15° frcm the stagnation point, the Reynolds number has 

to be increased to about 520 x 10^. The reference Reynolds number R c 

is based on the radius of the sphere j the boundary layer over the forward 
portion of a sphere is, therefore, highly stable at large Reynolds numbers 
with respect to the Tollmien-Schlichting type of waves. 


Because transition is often caused by roughness particles and because 
the roughness Reynolds number R^ where 


R h 



has been shown to be significant in determining whether a roughness par- 
ticle causes transition (see, for example, ref. 10), roughness Reynolds 
numbers were calculated. For roughness particles aman enough for the 
velocity in the boundary layer to vary linearly with distance in the 
region between the wall and the top of the particle, the roughness 
Reynolds number is 



However, since, from reference 3 or 4, 



then, 
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and because 0 = 


, I — , then 

V R c 


- -2 2 

D Ufh UhfR c 

\ = £• 

0v 


0 


R h = h 2 R c - 


3/2 Uf 

ft 


( 15 a) 


Because U, f , and Z are independent of h and R c , it follows that 


varies as h and as R 


3/2 


Equation (l5a) can also be written as 


K h = Uf 

h2R c 5/2 \[z 


(15b) 


Because the quantity Uf/VzT 
a plot of the quantity 


Rv 


for any value of R 


h 2 R c 
and h 


depends only on the pressure distribution 
enables R^ to be calculated rapidly 

for which the assumption of small height 




is valid. _ In_connection with the as sumption of smal l height it is noted 
that if h = 0 the error in u^ varies from about -20' percent for a 

velocity profile near separation to about l4 percent for the plane 
stagnation-point velocity profile; the error is zero for the flat-plate 
velocity profile. If, however, h = 215, then the error in varies 

from about -30 percent to about 30 percent for the same velocity profiles; 
the error is less than 2 percent for the flat-plate velocity profile. It 
is, thus, advisable that h not exceed 0 unless the pressure_gradient 
Is almost zero; in this case h can be slightly larger than 20. The 

Ry> _ . 

variation of — 7 — with x/c is shown in figure 4. For any fixed 

h 2 R c 5 / 

height h it i£ apparent that the maximum roughness Reynolds number occurs 
at a value of x/c of about 1 . 0 , or about 57 ° from the stagnation point. 

In order to give a clearer impression, figure 5 shows the dependence o-f R^ 
on x/c for a value of h of 0.005 inch on a sphere of 60 -inch radius at 
R c = 4.7 x 10°. 


Figures 1 to 4 allow the boundary-layer momentum thickness, stability, 
and roughness Reynolds number to be rapidly estimated for any sphere at any 
Reynolds number. These figures and figure 6, which gives the variation of 
the velocity-profile shape parameter K with x/c, can be used with the 
aid of references 3 and 4 to calculate all the other boundary-layer 
quantities . 



NACA TN 4350 


lit- 


Disk 

In order to investigate the effect on the various boundary-layer 
parameters of an increase in bluntness, the extreme case of the flow 
against a disk was investigated. Because no velocity distributions U(x) 
seem to be available for a disk, an approximate distribution was obtained 
by combining theoretical results with a measured distribution of U on 
the face of a plate of infinite span and perpendicular to the stream. The 
assumption was made that, at the same value of x. 



The theoretical expression for 


from reference 7 is 


(16) 


U; 


= - x 

* \Jl - X 2 


and the theoretical, expression for Up from reference 7 I s 


U 


x 

y/TT^ 


Thus, 



and 



U. 


p,a 


(IT) 


The values of Up^ a used were obtained from the experimental values 

given in figure 1 of reference 11. The variation of the approximate 

velocity distribution (5- ) , or U d , with x/c over the face of a disk 

VWd 

perpendicular to the stream is shown in f igure 7- A comparison of this 
velocity distribution with the calculated velocity distribution over the 
front face of the more forward of the two diBks that are separated by a 
cavity (see table III and p. 106 of ref. 12) indicates that the velocity 
distribution shown in figure 7 is probably a good-enough approximation 
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to the actual velocity distribution, to allow the main points to be shown 
by the boundary- layer analysis. In order to calculate the boundary-layer 
quantities of interest, a different and more rapid, although more approxi- 
mate, method was used to integrate the momentum equation for the disk. 

The method consists in obtaining an integrated form of the momentum equa- 
tion. The momentum equation (eq. (2)) can be written as 


M + £Z dr . gff . t - k ( 2 + . 0 

dx r dx U|_ 1 \ 0 / 


The quantity k is related to Z by the def ini tion 

k = Z^ 
dx 


From the relations 


and 


- 4 + 4 • *)] 


g 2 (K + 1) - fk x - k = 0 


(from ref. 4) It follows that 


f = 


. k ( 1 ~ I) 


+ £ 


2 ic 


z m - 1) 


+ g 


+ g 2 2 . 
S 6 


- f) 


+ g 


Now, let 




it . 
6 =b 


g - k^a = C (kj_ assumed to be constant) 


fr = d 


g^b _ 


( 18 ) 


(19) 


( 20 ) 


2 + V = A 
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where a, b, C, d, e, and A are constants. When the definition 
of k, the expression for f from equation (19) , and the definitions 
from equation (20) are substituted into equation (l8) , the result is 


dZ + ^ 1 dr + A - d dU 
dx r dx U dx 



( 21 ) 


In this equation 6*/ 0 , g, and k^ are assigned to be constant. Equa- 
tion (21) is linear in Z and is first order. Integration of equa- 
tion (21) gives 


Z = 


,ap2(A-a) 


a(e - kj J X ^ + 


2(A-d) 


( 22 ) 


When use is made of the fact that U = r n = 0 at the stagnation point 

p o u 

(x = 0) and that Z = 0 equation (22) becomes 


0 






(23) 


In order to find the value of 


approximate relations U 


- (Sk 


(e/iQt, note that near x = 0 the 
and r = x are valid. Equation (23) 


can then be integrated with the result that 




(24) 


Equation (24) indicates that, in the region in which the relations 
U = (g) x and r = x are approximately valid, the boundary-layer thick - 
ness is approximately constant. It is clear that in this region 0^R C is 
proportional to -7=^==. In figure 8 is shown the variation of — ^fT with 




x/c calculated by use of equations (23) and (24) together with k^ = 0 
and stagnation-point values for A - d and e; these values are 
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IT 


A 

e 


a = 3.3407] 


0.23002 j 


( 25 ) 


The variation of the boundary-layer Reynolds number Rq and the critical 
Reynolds number Rq c with x/c is shown, respectively, in figures 9 
and 10. The values of Rq c were calculated by first computing the 
values of k from 



p 

where Z = 0R C , then by reading the value of K for each value of k 

from a large-sized* plot of K against k, and then by reading the value 
of Rq from a large-sized plot of Rg c against K. Small plots of 

K against k and of Rg c against K are given in reference 4. 

The most striking results are that the boundary-layer thickness is 
greatest at the stagnation point and decreases with an Increase in dis- 
tance from this point and that the boundary layer is completely stable 

at a Reynolds number of 4.7 X 10^, based, on the radius of the disk. The 
smallest Reynolds number for which the boundary layer will first become 

unstable Is about 720 x 10^; the first unstable point will occur at about 
x = 0 . 75 • It is apparent that the boundary layer on a disk Is very sta- 
ble with respect to the Tollmien-Schlichting type of waves. 

In order to investigate the variation of roughness Reynolds number 

R h 

over the disk, the quantity = — j- was calculated. The results are 

p 5/2 

h R c 

shown in figure 11. In figure 12 is shown the dependence of R^ on x/c 
for a value of h of 0.005 inch on a disk of 60 -inch radius at 
R c = 4.7 x 10^. A comparison of figure 12 with figure 5 shows that the 
maximum value of R^ produced on a disk. is less than that produced on a 

sphere when the roughness height is the same for disk and sphere. A com- 
parison of figures 5 end 12 also Indicates that a roughness of given height 
produces a given roughness Reynolds number over a greater portion of the 
sphere than of the disk. This result follows from the fact that, for a 
given height of roughness and with the sphere and disk having the same diam- 
eter and Reynolds number, the maximum roughness Reynolds number on the 
sphere Is larger than that on the disk and occurs near x/c = 1.0 (near 57° 
from the stagnation point) on the sphere but at the edge of the disk. If 
the same roughness height is present on the sphere and disk and if the ref- 
erence Reynolds number is increased, a critical value of R^ will occur at 
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a lower reference Reynolds number on the sphere than on the disk and will 
occur farther forward on the sphere than on the disk. (See ref. 10 for 
information concerning critical values of R}*.) 


ACCURACY OF RESUETS 



In order to determine the accuracy of equation (23) , the quantity 
for the sphere was calculated for k, = 0 and with stagnation- 


point values of A - d and e; in figure 13 the results are compared with 
the results of the step-by-step integration of the momentum equation for 
the sphere. The agreement is very good not too near the separation point. 
For small values of x/c it is difficult to compute accurate values of 



by the use of equation (23) because both numerator and denominator 


become very small. In this region it is advisable to make use of equa- 
tion (24) as an aid in obtaining values of — |/r7. 

c ' c 


flow 


If the values of A - d and e 

( = 0 } are used, the values of 
\dx / 


for a flat plate parallel to the 
A - d and e, respectively, are then 


A - d = 3.5 


and 


e = 0.2145 


(ki - 0) 


With this value of A - d and the fact that r = sin x and U = ^ sin x, 

equation (23) can be integrated in closed form for the sphere with the 
result that 


= 


0.53492 


r 

105 


(sin x)^*5L584 


x - 


COS X 


8 


(sinTx + 


1 sirr^x + S s i n 3 x + 

6 24 48 


sin 


.1 1/2 


(26) 


A different closed-forma integral of equation (23), which is given in refer- 
ence 1, can be obtained from equation (23) by letting 
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a - a = 3 


and 


e = 0.2222 = o) 

The results of equation (26) are also shown in figure 13- Better agree- 
ment is obtained with the more exact step-by-step method of integration 
if the stagnation-point values of A - d and e are used in a region 
of falling pressure. Therefore, the values of equation (25) are used 
together with equation ( 23 ) to calculate the boundary- layer parameters 
on the disk. 

It may be worthwhile to note that equation (23) for the estimation 
of the boundary-layer momentum thickness is general in the sense that the 
constants e, A - d, and k^ are arbitrary. Therefore, the accuracy 

of a result can be improved by choosing values of the constants that are 
appropriate for the type of pressure distribution under consideration, 
(see table 2 of ref. 3 for suitable values.) One example of the effect 
of the choice of the constants is shown in figure 13- By making use of 
an average value of k^, equation ( 23 ) also permits a rough estimate to 

be made of the effect on the boundary-layer thickness of flow into or 
out of the surface. 


An indication of the accuracy of the calculated values of the 
boundary-layer thickness on both the sphere and the disk can be obtained 
from the exact solution for the boundary layer at the stagnation point 
of a body of revolution at zero angle of attack. This exact nondimen- 
sional velocity profile Is known to be the same as the velocity profile 

on the wedge with p = i (p. 129 of ref. 9). Following the notation of 

reference 2 and letting (3 - ^ gives 



the result Is 
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(27) 


The result obtained from substituting k^ = 0 and stagnation-point values 
from equation ( 25 ) into equation (24) is 



The error in equation ( 28 ) is about 7.5 percent. From this result it 
seems reasonable to infer that in the region of falling pressure the 
boundary-layer thickness, boundary-layer Reynolds number, and roughness 
Reynolds number are probably accurate to within approximately 10 percent. 

When, however, the accuracy of the critical Reynolds numb er Rq 

o , C 

is examined, the conclusion is reached that the calculated values of the 
critical Reynolds numbers can easily differ from the correct values by as 
much as 100 percent. Thus, from table II of reference 2 the value 2,662 
is obtained for R 0 ^ c at the stagnation point. The value of R 0 c 

obtained from the present analysis is 1 , 650 , a difference of about 60 per- 
cent. Both calculations of Rq^ c use the rapid method employed in refer- 
ence 15: the difference in the two values of R Q is caused by the dif- 

0 , c " 

ferent velocity profiles: one is the exact and the other is the approximate 

Schlichting profile. When it is noted that the rapid method of reference 13 
yields approximate answers, the conclusion is reached that the values of 
R 0 c given in the present analysis can easily differ from the correct 

values by as much as 100 percent. The conclusions reached in the analysis 
concerning the stability, however, still remain the same in a qualitative 
sense . 


CONCLUDING REMARKS 


The calculations for the sphere show that the boundary layer at the 
stagnation point of a sphere is much thicker than that on an airfoil, 
that the boundary- layer thickness increases very slowly with an increase 
in distance from the stagnation point, that the boundary layer over the 
forward portion of a sphere is highly stable at large Reynolds numbers 
with respect to the Tollmien-Schlichting type of waves, and that rough- 
ness of a given height produces the largest roughness Reynolds numbers 
at about 57 ° from the stagnation point. 
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The calc ula tions for the disk show the unusual result that the 
boundary-layer thickness Is greatest at the stagnation point and. decreases 
with an increase in distance from this point, that the boundary layer is 
extremely stable with respect to the Tollmien-Schlichting type of waves, 
and that roughness of a given height produces a given roughness Reynolds 
number over a smaller portion of the disk surface than over the sphere 
surface . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 9 , 1958. 
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Figure 2.- Variation of boundary-layer Reynolds number R 0 and of boundary-layer minimum criti 

cal Reynolds number FL „ with x/c for a sphere. 
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Figure 3 .- Value of Reynolds number R c 

at any value of x/c 


for neutral stability ^R Q 
on a sphere. 




Figure 5.- Dependence of Rh on x/c for £ = 0,005 inch at x/c on a sphere of 60-inch 

radius at R n = 4.7 x 106, 
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Figure 12.- Dependence of on x/c for h = 0.005 inch, at x/c on a disk of 60~ineh 

radius at R c 5= k.J x 10^. 




